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A scheme for addressing individual atoms in one- or two-dimensional optical lattices loaded with 
one atom per site is proposed. The scheme is based on position-dependent atomic population 
transfer induced by several standing-wave driving fields. This allows various operations important 
in quantum information processing, such as manipulation and measurement of any single atom, 
two-qubit operations between any pair of adjacent atoms, and patterned loading of the lattice with 
one atom per every nth site for arbitrary n. The proposed scheme is robust against considerable 
imperfections and actually within reach of current technology. 
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A system of cold neutral atoms trapped in an opti- 
cal lattice is expected to be a promising architecture for 
quantum information processing []]]. In particular, by us- 
ing superfluid to Mott insulator phase transition, a large 
number of atoms can be loaded into an optical lattice in 
such a way that each lattice site is filled with one atom 
Such a Mott insulator state naturally pro- 
vides a scalable qubit system, where each atom, i.e., each 
lattice site, stores a qubit encoded into its two ground 
hyperfine levels. Two-qubit operations in such a system 
can be achieved by exploiting cold controlled collisions 
between atoms @ or atomic dipole-dipole interaction 0] . 
Although these operations can not be applied selectively 
between one pair of atoms, they can be applied collec- 
tively between every adjacent atom, thereby generating 
a highly entangled state called a cluster state [8| . Once a 
cluster state is generated, arbitrary quantum logic oper- 
ations can be performed simply by measuring individual 
atoms in appropriate bases [9(. For this purpose, an es- 
sential requirement is the ability to manipulate and mea- 
sure single atoms individually. 

Addressing individual atoms in optical lattices is, how- 
ever, still experimentally challenging. The main reason 
is that the lattice separation of typical optical lattices is 
given by half the optical wavelength, while atoms are ma- 
nipulated by lasers having wavelengths comparable to or 
larger than the optical wavelength. This situation thus 
seems to imply that the addressing laser beam has to be 
focused within a length scale smaller than its wavelength, 
while it is limited by diffraction. 

There have been several approaches to this problem 
IIS 111 El- One straightfor- 
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ward approach is to increase the lattice separation, which 
was e mployed in most of the experimental attempts 



ISmii^lliS- In this case, however, it is hard to load a 
large number of atoms into the lattice with one atom per 
site. Another approach is to apply a position-dependent 
external field so that individual atoms are distinguished 
by their different resonance frequencies Hi EH- Oth- 
erwise, one would use a pointer atom that serves as a 
quantum information carrier In any case, the 

realization imposes significantly more demands on the 
experimental complexity and precision. 




FIG. 1: The involved atomic levels and transitions. Aj, 7i, 
and Qij (i = {m,q}, j = {1, 2}) denote the detunings, spon- 
taneous decay rates, and Rabi frequencies, respectively. 



In this paper, a new approach to the problem of ad- 
dressing individual atoms in optical lattices is proposed. 
One of the advantages of the present scheme is that it 
is compatible with any one- or two-dimensional optical 
lattice setup. For instance, it is compatible with a typ- 
ical lattice setup using an optical wavelength, and also 
compatible with a non-square (such as a triangular) op- 
tical lattice. Another advantage is that its experimental 
setup is quite simple: in addition to the original optical 
lattice setup, we need just several laser beams tuned to 
the atomic transition frequencies (with some detunings). 
Moreover, the scheme is robust against experimental im- 
perfections, and actually all techniques involved in it have 
been experimentally demonstrated. Remarkably, its ap- 
plicability is not restricted to the unitary rotation and 
measurement of a single atom. It can be also used for 
other important purposes, such as two-qubit operations 
between an arbitrarily chosen pair of adjacent atoms and 
the patterned loading of an optical lattice with one atom 
per every nth site for arbitrary n. 

Fig. Q] depicts the involved atomic levels and transi- 
tions. We use three ground levels |a), |6), and \q): among 
them, | a) and \b) are used to encode a qubit. Raman tran- 
sitions | a) <-> 1 6) and |6) <-> \q) are mediated by excited 
levels \e m ) and \e q ) with detunings A m and A q , respec- 
tively, where 74 and (i = {m, q}, j = {1,2}) denote 
the corresponding spontaneous decay rates and (com- 
plex) Rabi frequencies, respectively. The atomic state 
at lattice site s will be denoted by a subscript as | ) s . 

In order to illustrate the underlying idea, let us first 
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consider the case of a one dimensional optical lattice 
loaded with one atom per site. We consider the situa- 
tion where a single atom at lattice site k is manipulated 
by a well-focused laser beam centered on the atom. We 
assume that the beam waist is larger than the lattice sep- 
aration di , but is still smaller than 2Ldi for a moderately 
small integer L. The manipulating beam thus affects at 
most 2L + 1 atoms at sites k — L to k + L 18]. The basic 
idea is, prior to the single-atom manipulation, transfer- 
ring the atomic states of the neighboring sites coherently 
into quenching states that are decoupled from the manip- 
ulating beam. As will be shown later, this operation can 
be performed collectively for all neighboring 2L atoms 
simply by applying a couple of standing-wave fields. We 
will call this operation a quenching operation 24| . Once 
the quenching operation has been completed, the manip- 
ulating beam will affect only the selected atom at site 

(k) 

k. For convenience, let us define S L as the set of lat- 
tice sites s such that s = k + (L + l)n for all integers 
n. A quenching operation is defined such that it 
does not affect atoms at sites s £ while it transfers 
the atomic state coherently from the qubit subspace to 

the quenching subspace as \a) s (a\ + \q) s (b\ for all atoms 

( k) 

at sites s ^ S L , where state \q) serves as the quench- 
ing state. We also define the inverse quenching operation 
{Ql ) _1 transferring the state in the opposite direction. 

For the moment, let us assume we are able to perform 
the quenching operation. The applications of it follow 
rather straightforwardly. Some of them are listed below: 

Single-atom addressing. — The first one is a single- 



qubit unitary rotation R,i k \a) = exp(—iaa [ z h ' > /2), where 
0z fc ' = \a) k {a\ — \b) k {b\. In order to perform this op- 

(k) 

eration, we first perform the quenching operation Q L , 
then apply a focused laser pulse with Rabi frequency 
f2 m 2, and finally perform the inverse quenching opera- 
tion (Q\ If A m ^> |O m 2| 3> 7m, the laser pulse 
induces a unitary rotation R z k ^ ( \ n ™ 2 \ Atj, where At is 

the pulse duration. Note that this laser pulse has no 
effect on the neighboring atoms since after performing 
Qj^ there is no population in state \b) s for s $ . An- 
other single-qubit rotation about an orthogonal direction 
R { x k \(3) = exp(-i{3ai k) /2), whereat = |o) fc (6| + |6> fc (o|, 
can be performed with the help of a collective Hadamard 

operation He = (^) s °* , which can be done by 
applying global pulses addressing every atoms identi- 
cally. It is easily seen that a sequence of operations 
He, Rz (0), an d He results in the operation Ri'(f3). 
By combining these two kinds of rotations, an arbitrary 
single-qubit rotation can be performed. A single-qubit 
measurement can be also performed in the same fash- 
ion: it is done by observing the resonance fluorescence 
on transition \b) k <-» \e m ) k induced by a focused beam 
applied between performing and {Q^)~ l - 
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FIG. 2: Two plane waves incident on the optical lattice tilted 
respectively by angled 9n and 8i2 produce a standing-wave 
field along the trap axis. 



Two-qubit operation. — Another application is a col- 
lective CPHASE operation c£ fe) = <8> ae5 (»o (|o> a <o| <8> 

+ \b) s (b\ (gi ai s+1) ), where 1^ is the identity op- 
erator. This operation also starts from performing the 



quenching operation 



(fc) 



After that, we use the meth- 



ods outlined in Refs. [6|, [7| in such a way that the lattice 
potential corresponding to state |6) is shifted as much 
as about the lattice separation and restored back af- 
terwards so that every atoms in state \b) interact with 
the next neighboring atoms. Note that after perform- 
ing only atoms at sites s G S^' can be in state 
\b) s . This interaction can be used to perform a condi- 
tional operation such that the state acquires a phase ir 
only when it is |6} s \q) s+1 - By performing [Q^)~ l af- 
terwards, the net effect is the same as the desired oper- 
ation c[ fe) . If we use this collective CPHASE operation, 
even a selective two-qubit operation can be achieved. 
It is easily seen that a sequence of operations 



Rx k+1 \—ct/2), and C^ 1 results in a two-qubit controlled 
rotation |a) fc (ffl| <g> I {k+1) + \b) k (b\ ® R ( x k+1) (a) up to a 
single-qubit rotation R x k+1 \—a/2). 

Patterned loading. — Another straightforward applica- 
tion is the patterned loading of an optical lattice with one 
atom per every (L + l)th site. It is simply done by first 
preparing all atoms in state \b), then performing Q^ k \ 
and finally turning off the lattice potential correspond- 
ing to state \q). 

In the above scenario, the only nontrivial part is how 
to implement the quenching operation, which we will 
concentrate on in the remainder of this paper. For a 
quenching operation, we apply two standing-wave fields 
which, respectively, induce the transitions \b) <-> \e q ) and 
\q) <-> \e q ). Due to the position dependence of the fields, 
the transitions take place with different Rabi frequencies 
for each atom. Let us denote the corresponding Rabi 

(s) (s) 

frequencies at site s respectively by fl ql and fl q2 ■ The 

standing-wave field corresponding to £l q f (i = 1,2) is 
produced by two plane waves having a common wave- 
length \ qi and tilted from the trap axis by angles 9n 
and 6i2, respectively, as depicted in Fig. [51 If we set 
On = ±$i2 = 6qii the spatial period of the standing wave 



3 



along the trap axis is given by 



where the position 



of the nodes is determined by the relative phase between 
the two plane waves. What we should do is to adjust 
the involved parameters so that both the standing-wave 
fields have the same spatial period A = Xq i — Xq o 

* 1 H COSfql COS Uq2 

and have the common nodes. Furthermore, we require 

(k) 

that each lattice site in S L coincides with a node while 
all other lattice sites never coincide with the nodes. This 
requirement can be met in multiple ways, but for sim- 
plicity we restrict our consideration to the case of 



Xg/2 = (L + l)di. 



(1) 



To be more specific, let us assume the lattice poten- 
tial is sufficiently deep, thus the tunneling between lat- 
tice sites is suppressed. Furthermore, we assume every 
atoms are in the vibrational ground state. In this case, 
each potential well can be approximated by a harmonic 
potential. Let us denote by a s (a\) the annihilation (cre- 
ation) operator of the vibrational mode at site s. The 
involved evolution of the atomic state is then described, 
in the rotating frame, by the following Hamiltonian [l9j |: 

W« =(A q - i^) \e q ) s (e q \ + sm( V (a s + a\) + </> s ) 

x [Sli(t)\e q ),(b\+(h(t)\e q ),(q\ + h.c], (2) 
where r\ = 2Trxo/\ q is the Lamb-Dicke parameter with 



x (a s 



the position operator of the harmonic oscil- 



lation, <p s = 7r(s — k)/(L + 1), and Cli(t) (i — 1,2) is 
the corresponding Rabi frequency at the peak position of 
the standing-wave field. By using s (0| sm(r)(a s + a\) + 
4> s ) |0) s = e~ n / 2 sin(0 s ), where |0) s denotes the vibra- 
tional ground state of the atom at site s, the Hamilto- 
nian @ can be written as 



?#> H\-^)\e q ) s (e q \+e^ 2/2 



sin 



,) 



X [Sli(t) \e q ) s (b\ + n 2 (i) \e q ) s (q\ + h.c], (3) 
which implies the effective Rabi frequencies arc given by 



^ ) (t)=e-" 2 / 2 sin(0 s )O i (t). 
Consequently, we reach an essential condition 

fi# = for s e S { r k \ fi# ? for a & S ( L k) . 



(4) 



(5) 



(k) 

For L < 2, the quenching operation Q L can be per- 
formed straightforwardly by using the Raman transition. 



For a large detuning A q 3> 



n 



Raman transition 
)(^)o( s ) . 



\b) s <-> \q) s takes place with Rabi frequency Q ql 'fl q2 / A q . 
Since this Rabi frequency is the same for all lattice sites 



s $ when L < 2, the operation Q y £' is achieved sim- 
ply by applying a Raman 7r-pulse. However, this method 
does not work for L > 2 since in that case each atom 
experiences a different Rabi frequency. 



(fc) 




FIG. 3: Fidelity of state transfer |*} s = ^(|a) s + \b) s ) -> 
|*') s = ^(|a) s + \q) a ) by STIRAP. The fidelity between the 
final state and |$) s (|^') s ) i s plotted as a solid (dotted) curve 
with respect to the peak value fij^' of the Gaussian pulses. 



A more powerful way of performing the quenching op- 
eration regardless of L is to use the technique of stim- 
ulated Raman adiabatic passage (STIRAP) [20]. As a 
typical example of STIRAP, let us consider a case where 
the standing-wave fields are applied in order as Gaussian 
pulses 



(6) 



where the relevant parameters are chosen as t\ — t 2 = 6/j q 
and § t — 5/7 g . The detuning is chosen as A q = 1007 g . 
The resulting Rabi frequencies experienced by each atom 
have the same Gaussian pulse shapes but with a different 



peak value, which we denote by fl q — 



sin 



in accordance with Eq. ([J). A very useful property of 
STIRAP is that as long as Q q s ^ exceeds a certain thresh- 
old the state transfer is performed with a nearly ideal 
fidelity without introducing an unwanted phase change. 
In order to see this property, the fidelity of state trans- 
for I*)., = ^(|a) s + \b) s ) -> = ^(|a) s + |,),) by 

STIRAP has been calculated numerically with varying 
n q s) . In Fig. [31 the fidelities s (*|p s \^) s (solid curve) 
and s {^'\ Ps \^') s (dotted curve) are plotted, where p s is 
the final state after STIRAP. This figure clearly shows 
the desired result that the state \^) s is not altered when 

flq^ is small whereas it is transferred to |^') s when fl q s ^ 
is large enough. Consequently, the quenching operation 
can be performed with a very high fidelity if we 
take Qq so that fl q k+1>> is larger than the threshold (say, 
20-f q in Fig. [3]). Note that here we do not have to worry 

(s) (s) 

about the upper bound of £l q since for large fl q the 

(s) (s) 

dynamics is dominated by Q ql and fl q2 , in which STI- 
RAP depends only on the ratio between them, not on the 
individual amplitudes of them. 

Thanks to the property of STIRAP, the fidelity of the 
quenching operation is not sensitive to the exact shapes 
or amplitudes of the applied standing-wave pulses. How- 
ever, the fidelity would be affected more significantly by 
imprecise adjustments of the tilting angle 8 q i (i = 1,2) 
and the relative phase between two plane waves. Let us 
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denote the corresponding precisions respectively by A9 
and A(f>. Suppose that we require every lattice sites in 

(k) 

to coincide with the nodes within a spatial precision 
of Ad. Given a desired fidelity, the required precision Ac? 
can be determined from Fig. [3] For example, in order 

(k) 

for the fidelity at each lattice site s € S L to exceed the 
value corresponding to fl q = j q in Fig. [3] with choosing 



^(fc+i) _ 2Qj qi the required precision is given by 



Ad . 1 

V, 1 

d, ~ 20 k 



L + l 



(7) 



where sinc(a;) = 5HL£_ jf the total number of lattice sites 
is N, this precision is achieved when 



Ntaxi9 qi A9 + (L + 1 



A<j> 



< 



Ad 
~d7' 



(8) 



These inequalities give a trade-off relationship: in order 
to address a single atom in a larger optical lattice (larger 
N) with a less focused beam (larger L and 9 q i), we need 
more precise control of the tilting angle (AO) and the 
relative phase (A<fi) of the laser beams. For large L and 
\ q i /2 ~ di , we reach an asymptotic condition 



(L + l)(NA9 + A0/tt) < 1/20 



(9) 



For large N, this condition imposes relatively high preci- 
sion on the tilting angle. For example, if we take L = 4 
and N = 10 3 , we require A9 < 10 ^trad. This can be 
achieved through active stabilization based on a piezo- 
electic controller. 

Finally, several remarks are in order, (a) In case of 
two dimensional optical lattices, the present scheme can 
be applied in the same manner by introducing more 
standing-wave fields assigned to each trap axis, (b) 
The quenching operation starts from an assumption that 
there is initially no population in state \q). After one 



round of operation , manipulation, and 



(Ql 



(A\> 



however, some population unavoidably remains in state 
\q) due to the imperfections. In order to prevent this un- 
wanted population from accumulating, it has to be opti- 
cally pumped into the qubit subspace by using a globally 
addressing field after each round of operation (or pos- 
sibly after every n rounds for a moderate number n). 
Note that the optical pumping does not decrease the fi- 
delity of the preceding operation, (c) For an arbitrary 

single-qubit rotation, three pairs of and (Q^) -1 
are needed in general. It is because we use atoms having 
one quenching state as in Fig. [TJ If we introduce another 
quenching state \q') and perform the quenching operation 
such that Q ( £ ] = \q') s (a\ + \q) s (b\ for s £ S { L k) [El, an 
arbitrary single-qubit rotation can be performed in one 
round of operation Q ^ , manipulation, and (Q^) -1 . (d) 
A disadvantage of the present scheme is that all atoms 
have to be involved in addressing only one atom. For- 
tunately, the standard model of fault-tolerant quantum 
computation (FTQC) allows for such a situation since it 



regards both unitary operations and memories as being 
affected by the same amount of noise [22J . Together with 
global addressing, the present scheme is well suited for 
performing collective operations of the form (^) seS (fc) 

and <g) seS (fc) V^ s ' s+1 \ where and V^ s ' s+1 ^ are one- 
and two-qubit operations acting on lattice sites s and 
s + 1. If we take L + 1 as the number of qubits in one 
error-correcting block, such collective operations are ac- 
tually what should be done in the error-correcting steps 
of FTQC. In this context, further investigations on the 
effects of reduced parallelism and restricted geometry 
would be desirable. 
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